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Abstract 



qq ■ This paper examines a heterogeneous beliefs model in which there is a process 

CN \ that is only partially observed by the agents. The economy contains a risky as- 

set producing dividends continuously in time. The dividends are observed by the 
agents. The dividends are assumed to be a known function of some other unob- 
served process. The agents use filtering to estimate the value of this unobserved 
process. The agents have different beliefs about the dynamics of the unobserved 
process and therefore form different estimates. We analyse this model and derive 
the state price density. We use this to derive the riskless rate. We also charac- 
terise the price of the risky asset in terms of the solution of a series of differential 

equations, 
i — i ■ 

^> 

in ! 1 Introduction 

On 



This paper will build on the theory outlined in iBrown &: Roger 1 (|2009ah . We 



consider an economy with unit supply of a risky asset and zero net supply of a 
riskless asset. The risky asset produces dividends continuously in time and these 
dividends are observed by all the agents. The dividend process is assumed to be 
a function of another process, X, which is not observed by the agents. In order 
for the agents to determine how to behave, they will use filtering to estimate the 
$_i ' unobserved process. Since the agents never observe this process, we will obtain a 



non-trivial steady state. 

The model assumes that the dividend process is a linear function of the first 
component of the process X; this linear function is known to all the agents. We will 
assume that X is a multidimensional Ornstein-Uhlenbeck (OU) process. Since X is 
a multidimensional OU process, all the components of X will potentially affect the 
behaviour of the first component. Agents can easily work out the first component 
of X by looking at the dividend process, but they will be unable to work out the 
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other components. Since they need to know the other components to predict the 
behaviour of the dividend process, they will need to use filtering. 

We will assume that the agents have different beliefs about the stochastic differ- 
ential equation that governs X. Specifically, the agents each have a different belief 
for the matrix that determines the dynamics of the OU process. They will therefore 
each compute different estimates of X. This will therefore affect how they behave. 

The agents maximise the discounted expected utility of future consumption. 
They each have CA RA utility. We can use the optimisation problem to characterise 
the equilibrium as in Brown Sz Rogers] ( 2009a ). Using market clearing we may then 
determine the state price density. We then use this to work out the interest rate 
process. We can also explore the stock price, although we are unable to obtain a 
closed form expression. 

Ther e is a large literature on heterogeneous b eliefs, which h a s been discu s sed in 
detail i nlBrown fc Rogers) (l2009ah . Work includes iKurzl (l2008bh . iKur d (Jl994h iKurzl 



(1997), Kurz & Motolese (2006), Kurz (2008a), Kurz et al. (2005), Fan 



Harrison Krepd (Il978h. iMorrisI fll996lV. |Wu fe Gucl (l2003ll. Iwu & Guc (2004), 



Harris Ravivl (|l993l ) , iKandel Sz Pearson! (|l995l ) , iBuraschi Sz Jiltsovl (|2006l ) , iJouini Sz Napp 
(120071 1. We will now briefly review some of the filtering models that are of a similar 
fla vour to the on e prese nted here. 

Xiong &z Yan (|2ooah consider a model in which there is a risky technology whose 



instantaneous return, f t , obeys an OU process, but the mean of this OU process is 
not fixed. This mean is unobserved by the agents and follows a different OU process. 
Agents observe / and another information process which is in fact independent of 
all the other processes in the system. However, agents think the information process 
is correlated to the Brownian motion driving the unobserved OU process. The two 
groups of agents have different beliefs about whether this correlation is positive or 
negative. The paper shows that trading volume increases with belief diversity. The 
paper also looks at bond yields and use s the m odel to explain some effects that 
are observed in bond data. Dumas et al. (|2005h consider a similar model in which 
agents observe a dividend process which is obeying some SDE. The agents do not 
know the drift of this SDE. This drift is random and obeys some other SDE. Agents 
also observe an information signal. One group of agents thinks that this signal is 
positively correlated with the dividend process. The other group correctly thinks 
that it is not. Agents must perform filtering to deal with the unknown drift. The 
authors derive asset prices and show that the presence of irrational agents leads to 
lo wer stock prices that ar e more volatile. 



Brennan &z Xial (120011 ) consider a representative agent model in which the div- 



idend process of a stock satisfies an SDE that has an unobserved drift. The drift 
satisfies an OU process. The agent maximises the expected discounted utility of 
future consumption. The authors solve for the stock price and hence deduce quan- 
tities such as the equity premium and price dividend ratio. They then calibrate the 
driving SDEs of their model to the Shiller data set and calculate various derived 
quantities. Their results provide good agreement with the data, but they assume a 
constant of relative risk aversion of 15 and negative discounting, both of which are 

un desirable. 

Scheinkman k Xiond (|2003h consider a model in which two groups of agents 



observe a dividend process that obeys an SDE with drift unknown to the agents. 
This drift obeys some other SDE. The agents observe two signals, called A and 
B. Agents in group A think that the noise term in the SDE driving signal A is 
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correlated to the noise term in the SDE driving the drift process. Similarly, agents 
in group B think that the noise term in the SDE driving signal B is correlated to 
the noise term in the SDE driving the drift process. In fact, the noise terms in 
the SDEs are independent. Using standard filtering theory, the authors deduce the 
conditional means of each of the agents for the drift term. They interpret this as 
the beliefs of the agents. The authors derive equilibrium prices for the asset. They 
then show that each agent is willing to pay more than their perceived fundamental 
value of the asset, since when they buy the asset, they are also buying an option to 
res ell the a sset a t a possibly inflated price. The authors use this to explain bubbles. 
Davidl ( 20081 ) introduces a model in which there are two fundamental processes 



in the economy. The first is the output of the economy, of which each agent receives 
a fraction. The second is a dividend process. The drift of the dividend process and 
the economic output obey two SDEs. The drifts of these SDEs are unobserved 
by the agents and follow a two-state continuous time Markov chain. Agents form 
beliefs about the drifts based on their observations. They then seek to maximise 
their expected utility. David derives quantities such as the consumption and in- 
terest rate and uses his model to explain why the stock volatility can be so high. 
David Veronesil (2002) also look at a continuous time model in which at any given 



time, the economy can be in one of two states; boom and recession. The agents do 
not observe this state directly, but instead must infer it from their observations of 
the dividend process. 

Detemple fc Murthvl (| 19941 ) consider a model with 3 different processes. One 



process gives the total amount produced by the economy. This process is observed. 
However, its drift is unobserved and obeys a different SDE. Finally, there is a third 
information process, observed by the agents. The drift of this process is unobserved. 
Agents work under different measures and form beliefs about the unknown drifts. 
Agents maximise expected utility. The paper derives an equilibrium for this model. 

The structure of this paper is as follows. Section [2] explains the model in detail. 
Section then proceeds to analyse the model by deriving a state price density. In 
section H] we look at the riskless rate and stock price. We can proceed far with the 
calculation of the stock price, but unfortunately are unable to obtain a closed form 
expression. Finally, we conclude in section [5] 



2 The Model 



We will now describe our model in detail. As in Brown &; Roger 1 (|2009ah . 



we 

assume that there is unit supply of a single risky asset. This asset produces a 
dividend continuously in time. There is also a riskless asset in zero net supply. The 
dividend of the risky asset at time t is denoted St- We will assume that 5t is a 
function of the unobserved process Xf. Agents disagree about the dynamics of X. 
Thus, when they perform filtering they will form different estimates of X, which in 
turn will affect their behaviour. 

We take our sample space to be = C(R + ,M n ) and we let X t (uj) = u(t) denote 
the canonical process. Let Po denote the measure under which A is a standard 
Brownian motion in M n . 

There are J > 2 agents in our model. Each agent will work under his own 
measure. Agent j thinks that the true measure is given by P J , where P- 7 is defined 
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by: 



A? := - 



exp 



(B 3 X s ,dX s )-- / \B 3 X s \ 2 ds 
o * Jo 



where B J is an NxN matris0. The Cameron-Martin-Girsano^] theorem gives that: 

dWl :=dX t + B j X t dt 

is a standard Brownian motion under P J . Thus, X is an OU process under ¥ J . 

None of the agents can observe X directly, but they do observe the dividend 
process, St- The dividend process is given by: 

6t = a + cei • X t 

for some constants ao and a, which are known to the agents^]. We assume that all 
the agents receive the same signal, namely the dividend process, St- The important 
point is that in general B l ^ B 3 for % ^ j, so their beliefs will differ. For example, 
agent 1 may take B 1 to be diagonal and hence 5t will behave as a one-dimensional 
OU process. In contrast, agent 2 may construct B 2 in a way that means 5t is 
affected by all the different components of Xt. There are many different possible 
setups, but the exact form of the beliefs is not important at this stage. 

We define the observation filtration to be 3^t := cr{e\ -X s : s < f ) . The objective 
of the agent j is: 



max 



Uj(t,4)dt 



where cj- is the consumption of agent j at time t and Uj(t,cj.) = ——e 7jX e pt is 
the discounted CARA utility. Alternatively, we may write the objective as: 



max E 1 



(2.1) 



where K° t is defined by: 



^:=E [A^|3V 



(2.2) 



This completely specifies our model and we can now proceed to deduce a state price 
density. 



3 Deriving the State Price Density 



D.T 

paper of iBrown Roeersl (|2009bT ). Thus, we may proceed in exac tly the same way 
to deduce the state price density. We do this in two steps (see iBrown &; Rogers 



1 We will assume that all the eigenvalues of B 3 have neg ative real part. This implies t hat under the 
measur e P J , X is stable in distrib ution (see Corollary 2.2 of Basak fc Bhattacharva (1993)). 
2 See lRoeers fc Williamsl (|2000l ). IV.38 for an account 
3 Here, e\ is the vector with 1 in the first position and zeros elsewhere 
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(|2009al ) for more details). Firstly, by looking at the price of an arbitrary contingent 
claim we can deduce that: 



where i>i is some some constant, different for each agent. Taking logarithms and 
using market clearing, we may then deduce that: 

logCt = G — pt — TaX t + 7 ^ - log Af 

where T _1 = YljlJ 1 an d G is some constant. We can use this state price density 
to deduce the interest rate and the stock price. 

3.1 Calculating A t 

In order to understand the state price density, we need to obtain an expression for 
A-j. For the analysis that follows, we will omit the superscript j unless doing so 
would cause confusion. 

First note that 3^ is generated by the process X 1 and furthermore that X 1 is a 
3^-Brownian motion under Pq. Since At is a ^-martingale, we may therefore apply 
the martingale integral representation theorem^ to deduce that: 

A t = 1 + f H s dXl (3.1) 



o 



for some ^ process H t . 

We will now determine the process H. To do this, we introduce an arbitrary 
process 9t, which is assumed to be bounded and adapted to the filtration (3 ; t)t>o- 
We have: 

E°[At f 9 u dX l u ] =E°[At f e u dxi} 
Jo Jo 

= K°[[ A s (-BX s ,dX s ) [ 9 u dXl) 
Jo Jo 

= E°[[ K,e u e' x (-BX u )du] 
Jo 

= E°[f e u e[(-KBX u )du] 
Jo 

where Y u denotes E°[5 / U |3^n] and e% denotes the column- vector with 1 in the i- 
th position and zeroes elsewhere. Here, we have used the definition of conditional 
expectation in the first and final lines, the integral representation of A in the second 
line and the fact that E.AB = EL4, B] in the third line. However, we may also use 
(|3.1|) to deduce that the above expression is: 



E°[/ H u 9 u du] 
Jo 



See lRogers fc Williams (2000) IV.36 
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Thus we obtain: 



But since 



H)t>o 



is just an arbitrary (3^t)t>o-pi'ocess, we deduce that 
H u = -e[(KBX u ) 



E°[A„|^] 
—e'iBX u A u 



where E denotes expectation under the agent's measure. Combining this with (|3,ip . 
we have: 



dA u = -A u e[BX u dXl 
We therefore need to develop X in order to better understand A. 

3.2 Deriving X u 

First define: 



(3.2) 



Then we have that: 



Also define: 



h\ = -ei(BX t ) 



dX\ = dW} + h\dt 



N t = X\ — I h\ds 



(3.3) 



Furthermore, define the conditional covariance matrix V by: 

v? =xjxi- x]xi 

We wish to derive the distribution of Xt- Note that since [Xl,Xt) is a Gaus- 
sian process, the conditional law of X t given y t will be Gaussian. Thus, we just 
need to determine the conditional mean and conditional covariance matrix. It is 
straightforward to deduct that: 



X\ = Xq — {BX s )ids - / ((BV)u - 5n) dN t 



(3.4) 



Alternatively, in vector notation: 

X t = X - f(BX s )ds - I (((W)').i - e x ) 
Jo Jo 



dN, 



3 See the appendix 
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We may also deduce that: 



dV s ij = - [(BV)ji + (BV)ij - Sij + {{BV) U - 5 a )((BV) lj - 5 n )} ds (3.5) 



Note that there is no dN term in the above SDE and thus V is in fact deterministic. 
We may also use (|3.5p to specialise to the case in which V = 0, but we first introduce 
some notation. 

3.2.1 Notation 



Since = V^ 1 = 0, we may write: 



Vt 



0„-l V t 



where ra _i denotes the zero- vector of length n — 1. Similarly, we may write: 

B : 



hi c 
A B 

Furthermore, we may define xt and xt by Xt = (xt,xt)'. We are of course only 
interested in working out xt- From (13. 4D we deduce that: 

dx t = -{Ax t + Bx t )dt - (VC)dN t 

But, 

dN t = dx t + Qrnxt + C'x t )dt 

so we learn that: 

dx t = -(A + VCb u )x t dt - VCdx t - (B + VCC')x t dt (3.6) 
3.2.2 Stationarity 

We now can examine the steady state in which V is not changing. We want to solve 
V = 0, or equivalently: 

[(BV)ji + (BV)ij - 5^ + ({BV) U - 8ii){{BV)x j - 8 jX )] = (3.7) 

This equation is automatically satisfied whenever i = 1 or j = 1. For i, j > 1, we 
have: 



(BV)ij = ^BikVkj = ^J-Bi(jfc-i)%_i)j = BVi 

and: 

BV ll = B ikV kl = C k V ki = (C'V)i 

k k=l 

Substituting in (|3.7p we see that the stationary covariance matrix is specified by: 



k=X k=2 



n-X 



BV + V'B' + V'CC'V = J, 



n-l 
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3.3 The State Price Density 
3.3.1 Notation for multiple agents 

Now that we have developed our expressions for A and X, we may use these to 
shed light on the state price density. 

Before we do this, we define some further notation to deal with the J different 
agents. Recall that each agent has a different B matrix, and hence a different B, C 
and A matrix. Let x\ be the estimate of (X^,Xf,...,X^)' according to agent j. 
Then we may define: 

%t = i&ti^ii &t) 

to be the concatenation of the estimates of all the different agents of the unknown 
components of X. Using equation (13. 6|) we see that Zt satisfies the SDE: 

dZ t = (Aix t + B\Z t )dt + Q\dxt 

Here, A\ is the stacked version of —(A + VCbu) for each agent and Q\ is a stacked 
version of —VC for all the different agents (recall that these matrices are different 
for each agent). B\ is the ((n — l)J)x((n — 1) J) block diagonal matrix with each 
of the —(B + VCC') of each of the different agents in the J blocks. 
We also define Zt = (xt, Zt)' . Then we have that: 



dZ t = BZ t dt + Qdxt 



where Q = (1, Q^)' and 



B 



Oj(n-l) 



3.3.2 Final expression for the state price density 

Having defined this notation, we may now use it to derive our final expression for 
the state price density. 

We know from (^]) that the A^ satisfies the SDE: 

dA{ = -Kie'^XidXl 
The solution to this SDE is given by: 



Ai = A J n 



b ex p 



[ -{^)Xidx u -\ ! (e'^Xifdu 
Jo z Jo 



In view of the expression we have for the state price density, we will need: 

— log A J T = — log A J + / -Xidx u - - / ^ ^du 

7j 7i -A) Ij 2 Jo 7j 

= — log A J + / ( - -^-x u - ± — L xi)dx u - - / 7i(— ^ + - — -a&) du 



Specifically, we are interested in the sum of these terms, which can be expressed as: 



V — log A^ = const. + [ ( - V ^-x u - V 

f-T 



7j 



m r 2 



b> n (C>)' 



,1 .,/.?.' 



= const. + A a Z t dx t — — f Z t 'j3Z t dt 
Jo 2 Jq 

where a is a 1 + J(n — 1) dimensional vector given by: 



, +E (4)<cy)'S„ 



7j 



7j 7i 



72 



7J 



and 



/y- Mil! fe n( g1 )' b n( g2 )' 
7j 71 72 

^ -yj u n— l,n— 1 







6^ 

72 



On- 



l,n-l 



C 2 (C 2 )' 



72 



V 



"ll^ r\ n 

j u rt— l,n— 1 u n— l,n— 1 



7J 

Un— l,n— 1 
O n -l,n-l 



C J (C J )' . 
7J / 



where O n _i >n _i denotes the zero matrix of dimension (n— l)x(n— 1). Note that in the 
above notation the first row and column have dimension 1, whereas all subsequent 
rows and columns have dimension n — 1. 

The final expression we obtain for the state price density is given by: 



log ( T = const. - pT- Tox T + - ( / a'^cixj - - / Z t 'pZ t dt) 

J Jo 2 Jo 



4 Using the State Price Density 

We now have a very explicit form for the state price density, which can be used to 
price any assets. We will derive the riskless rate in this model and also illustrate 
the method for calculating the stock price. 



4.1 Riskless rate 

Recall that there is a single riskless asset in the model, which is in zero net supply. 
The interest rate must satisfy: 

d(t = (t(-ndt - Ktdx t ) 

where (r t ) t >o is the interest rate process and (K t )t>o is some other process, not 
currently of interest to us. Applying Ito's formula to our expression for the state 
price density, we may simply read off the interest rate process as: 

r t = p + ^jZ^Z t -^(ja'Z t -a) 2 
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4.1.1 Remarks on the riskless rate 

We see that if the agents are impatient (p large) then the riskless rate offered to 
them must be higher to stop them simply consuming their wealth. If we consider 

r 

the case in which Z t = 0, we see that the riskless rate is simply p — —a" ! , so that a 

larger volatility in the dividend process means that the riskless rate is lower. The 
dependence on Z is more complex, since it will in turn depend on our assumptions 
for a and (3. 

4.2 The Stock Price 

The stock price is given by: 

s t =K[r^dT] 
jt « 

/oo 
e - P (T-t) E o [ exp ( _ Ya{x T - x t ) 

r f T i r T 

— {J a'Z u dx u --J Zu f3Z u du))(ox T + a )]dT 



+ 

If we can calculate: 



E? 



r - i 

[exp ( - T(j(xt - x t ) + — ( J a'Z u dx u --J Z u '(3Z u du) + 6(ax T + a ))] 

then we may differentiate with respect to 6 to obtain the conditional expectation 
that we require. 

4.2.1 Calculating the conditional expectation 

We will calculate: 

V T (t,Z t ;9) :=E t [exp{-ra(x T -x t ) + j(J a'Z u dx u 

Z u 'pZ u du) + 9(ax T + oo)}] (4.1) 

We will show that: 

V T (t, Z t - 9) = exp (±(Z t )'a(T)Z t + b(r)Z t + c(r)) 

where r = T — t and a(r) is a symmetric ( J(n — 1) + l)x( J(n — 1) + 1) matrix, 
6(r) is a (J(n — 1) + 1) row vector and c(r) is a scalar. We will omit the explicit 
dependence on r except where we specifically require it. In order to calculate a, b 
and c, we will use a martingale argument. We define: 

M? = E°[exp{-IVx T + -J f a'Z u dx u -\ f Z u 'pZ u du) + 9(ax T + o )}] 

J Jo 1 Jo 

= V T (t, Z t ; 9) exp [ - Tax t + %{ I a' Z u dx u - \ [ Z u 'pZ u du) 
1 J Jo 1 Jo 
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Since (M^)o<t<T is a martingale we will apply Ito's formula to it and deduce that 
the term in dt is zero. This will enable us to get a series of three differential 
equations for a, b and c. Firstly, we have that: 

d(^Z[aZ t + bZ t + c) = ( - -Z[aZ t -bZ t -c + bBZ t + Z[aBZ t + -Q'aQ)dt 

+ (bQ + Z'aQ)dx t 

Thus, 

^ = ( - ^dZ t - 6Zt - c + 65^ + Z t 'a5Z t + ^Q'aQ + ~(6Q + Z'aQ) 2 )dt 

+ (bQ + Z'aQ)dx t 

Now define: 

K(t,Z t ) = -Tax t + — ( / a'Z u dx u --l Z' u f3Z u du) 
J Jo £ Jo 

Then, we have 

d(exp{i^(t,Z t )}) r /«% \ t , ( r =,5= 1 r2 ,a'Zt \2\ 

- r(— — - <r)da; t + ( - —:Z t 8Z t + -T (— — - a) jrft 



exp{K(t,Z f )} J V 2J 2 v J 

So finally we may calculate cfM: 

= ( - ^ZjdZt - 6Z 4 - c + + + ^Q'aQ + i(6Q + Z'aQf)dt 

+ ( " 27^* + ^ r2 (^r " rf) dt + r( ^r " + 2 ' a Q) dt + i-> dx * 

Since M is a martingale, the dt term in the above equation must be zero. Thus we 
obtain: 

1 - T - -- T 2 2TaOa' 
-Z[{-a + 2aB - -3 + aQQ'a + -=aa / + 

2 J J z J 

+ {_6 + 6 £ + 6QQ' a - 1^ + ™^ - rVQ'a}^ 

+ {-c + \Q'aQ + ^QQV + ^rV - IV6Q} = (4.2) 

Each of the above expressions in the parentheses must be equal to zero. We also 
have some boundary conditions; note from (|4.ip that: 

V T {T, Z T ; 9) = exp (9ax T + 9a ) 

Thus, we have the boundary conditions: 

a(0) = Oj(„_i),j (n _i) b(0) = (9a, j(b _i))' c(0) = 0a o 

Unfortunately, the first line of (|4.2[) gives us a matrix Riccati equation. It seems 
that solving this equation is intractable, and thus we are unable to proceed further 
with the stock price. 
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4.2.2 A PDE approach 

An alternative method for tackling the stock price is to use a PDE approach. We 
may argue that QSt + / * ( s o~ s ds is a martingale. We can then apply Ito's formula 
and deduce that the dt term is zero. This will give us a PDE for the stock price. 

Proceeding in this way, first define St = h(Z t ). Then we may deduce that the 
stock price must satisfy the PDE: 

Vh(Z t ) ■ (BZ) + ^Q'H(Z t )Q + {ja'Z - Ta)vh{Z t ) ■ Q - r t h(Z t ) + 5 t = 

where H denotes the Hessian of h. Unfortunately, there seems to be little hope of 
solving this PDE, given that the dimension of this problem is J(n — 1) + 1. 

5 Conclusions 

We have introduced a very general model for dealing with heterogeneous beliefs of 
agents. This model assumes that there is some unobserved process X that drives the 
dividend process. Agents differ in their views about the SDE that X obeys. This 
affects their behaviour. We are able to proceed far with the analysis of this model, 
in particular, deriving the riskless rate in this model. However, the calculation of 
the stock price appears intractable because it requires a solution of a matrix Riccati 
equation. 
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APPENDIX 



A Deriving expressions (13.41) and (13.51) 



We seek to determine an SDE for X . In order to do this, we will follow section 
VI.8-9 of Sogers fc Williamsl (|200Ch . 

First recall that N was defined in by: 



N t = X\ 



h\ds 



By Theorem VI. 8.4 of lRogers &: Williams l|200d ) we know that Nt is a 3^ Brownian 
motion. 

Let / : lZ n i— > 7?. and let /t := f{X t ). Let (7,/t be the function such that: 



M t = ft- /< 



•-/ 



Gfsds 



is an J^-martingale. Further, let aj be such that: 

dftdXj = atdt 

Then by (VI. 8. 16) of lRogers fc Williamsj l|200d ). we have: 



ft = fo+ f Gfsds + / ( /.//J - Mi + nl ) dN s 
Jo Jo 



(A.l) 



We will use (]A. 1[) to derive the distribution of A t . Note that since {X},X t ) is a 
Gaussian process, the conditional law of X t given y t will be Gaussian. Thus, we 
just need to determine the conditional mean and conditional covariance matrix. We 
do this by applying (jA. 1[) to the functions: 

• nf | v -v* ■ T* ■ 

. JU 1 7 it' 'l it' j 

Firstly, working on the /*, we have: 

dX\ = dWt + h\dt 
Noting that df\dX\ = dW\dWl = 6 a dt, we have that: 

Qft =hl a l = 5 a 
If we put all this in (jA.ip . we obtain: 

Xl = 5? -J e[BX s ds - J e' x (aJbA, - X\J x BX a - 5 tl ) dN s 



Upon substituting for V, we obtain (|3.4p . 
We now move onto f lJ . First note that: 



d(xjxi) = xidwf + xldwi + [X*h> + X>h* + <5y]dt 
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and that: 

dffdXl = XlSjxdi + XfS a dt 

Thus, we obtain: 

Gff = Xihl + x{h\ + % a? = XtSj! + X'Six 
Applying (jA.ljl . we obtain: 



X\X\ = X^Xl + J \X\b? a + Xihl + 4, ) ds 

ft 



+ 







XjXihl - X\Xih\ + XjSji + Xl5a 



dN, 



x o x o 



(BX s )jXi + (BX s )iXl - 6ij ds 



E B lk X*XiX j s - E B Xk X*Xixi - X\5 jX - XiSn 



dN s 



We need to work out terms of the form X^XiXg . However, since X is a Gaussian 
process, we have that: 

x*xix* = xjv s jk + x j s vi k + x}v s ij + xjxix} 



If we now define := X\X\, we obtain: 



x\x\ = x^xi - / E B i* x s x i + E B * x i x i - 8*3 ds 

J ° \ k k / 

- /* [E B ikXiv. ik + E B ^nv: k + E b^xw 

J ° k k k 

+ E B lk XiX J s x} - E B lk x}xixl - Xi5 3l - xi& 



x o x o 



! ( E b m v + M )ki + E B *<y + M )u - s a ) ds 

- j [(BVhjXi + (BV)uXl + {BXs^V? 

- (BX^V^ -Xidj! - Xid a 



Thus, 



dXiXi = - [{B{V + M))ji + (B{V + M))ij - Stj] ds 

- [{BV) Xj Xi + {BV)uX j s - X\S jX - XiSa] dN s (A.2) 
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We now proceed to calculate: 

dV l s j = d(XiXi - X]xi) 

Noting that 

dXl = -{BX s )ids - ({BV)u - 5 a )dN. 

we see that: 

d(xlX J s ) = [ - Xt(BX s )j - Xi(BX s )i + ((BV)u - 5 ll )((BV) lj - S^ds 

- [Xi{{BV) Xj - Sjx) +xi({BV) li - 5 n )]dN s 
Putting this together with (|A.2|) gives: 

dV s ij = - [(B(V + M))ji + (B{V + M))ij - Sij - Xi(BX s ) 3 - xi(BX s )i 

+ ((BV) li -S a )({BV) lj -5 jl )]ds 

But, we note that: 

Xi{BX B )j =Y J B jk xlxi = {BM)ji 

k 

so finally, we deduce (|3.5p as required. 
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